MOCK GAUSSIAN BEHAVIOUR FOR LINEAR STATISTICS OF 
CLASSICAL COMPACT GROUPS 



CP. HUGHES AND Z. RUDNICK 



Abstract. We consider the scaling limit of linear statistics for eigenphases 
of a matrix taken from one of the classical compact groups. We compute 
their moments and find that the first few moments are Gaussian, whereas the 
limiting distribution is not. The precise number of Gaussian moments depends 
upon the particular statistic considered. 



1. Introduction 

In this paper we investigate the scaling limit of linear statistics for eigenphases of 
matrices in the classical groups. Given a unitary N x N matrix U with eigenvalues 
e l8 '\ 1 < n < N, and a test function g which we assume is 27r-periodic, consider 
the linear statistic 

N 

Tr g(U) :=X>(0«) 

ra=l 

A number of authors have studied the limiting distribution as N — > oo of Tr g(U) 
as U varies over a family G(N) of classical groups and have concluded that the 
distribution is Gaussian, see |2j [j], f|. 

Soshnikov || showed that this result remains valid in the "mesoscopic" regime, 
that is if one considers eigenphases 9 n in an interval of length about 1/L where 
L = Ln — » oo but L/N — > 0: For a Schwartz function / on the real line, define 

OO j 

Fl(0):= f(^(0 + 2nj)) 

j=-oc 

which is 27r-periodic and localised on a scale of 1/L. Soshnikov || showed that as 
long as L/N — > 0, then the limiting distribution of Tr Fl(U) as U ranges over all 
unitary matrices in U(iV), N — > oo is a Gaussian with mean 



N f°° 

L J J {X) dx 



and variance 

f(t) 2 \t\dt 
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where the Fourier transform is defined as 

/oo 
f(x)e~ 2 ^ t dx 
-oo 

There are similar formulae for the other classical groups. 

Our goal is to investigate these linear statistics in the scaling limit, that is to take 
L = N. Thus we set 

N 



Z f (U) := Tr F N (U) = £ F N (9 n ) 



In P| we proved 



Theorem 1. If supp / C [— 2/m, 2/m] then the first m moments of Zt{U) over 
the unitary group U(iV) converge as N — > oo to the Gaussian moments with mean 
I— f{ x )^ x an d variance 

/•OO 

min(M,l)|/( U )| 2 d U 

We called this a "mock-Gaussian" behaviour. It is worth remarking that in S we 
find the full distribution of Zf, and it is not Gaussian. Only the first few moments 
are. 

The purpose of this paper is to demonstrate mock-Gaussian behaviour for linear 
statistics in other classical compact groups, the special orthogonal group SO (TV) 
and the symplectic group Sp(iV) (N must be even in the symplectic group). If 
e 10 is an eigenvalue of a matrix U taken from one of these groups then e~ l6 is an 
eigenvalue too. This means 1 is always an eigenvalue of U 6 SO(N) if TV is odd. 

Due to the pairing of eigenvalues, the function / must be even. Our results are 

Theorem 2. i) //supp/ C [— 1/m, 1/m] then the first m moments of Zf(U) over 
the symplectic group Sp(A^) converge to the Gaussian moments with mean 

7(0) - / /(«) du 
Jo 

and variance 

2 / \u\f{u) 2 du 

J -1/2 

ii) //supp / C (—1/m, 1/m) then the first m moments of Zf(U) over the special 
orthogonal group U G SO(N) converge to the Gaussian moments with mean 

/(0) + f f(u) du 
Jo 

and variance 

2 / \u\f{u) 2 du 

J -1/2 

Remark. There exists / such that supp/ C [—1/m, 1/m] and whose m + 1-st 
moment is not Gaussian. 
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1.1. Moments and cumulants. One approach to proving such results is to use 
the Fourier expansion g{9) = ^2 n g n G and expand Tr g(U) as a sum 

Tr 5 ([/) = ]T.g„Tr([/") 

n 

Computing moments of Tr g(U) then boils down to being able to compute integrals 
of products of Tr([/ n ) over the classical group. Theorem|l]for the unitary group was 
proven in || using this approach by employing a result of Diaconis and Shahshahani, 
|2j ||, concerning moments of traces of random unitary matrices. Their result 
is a consequence of Schur duality for representations of the unitary group and 
the symmetric group, and the second orthogonality relation for characters of the 
symmetric group. 

The paper by Diaconis and Evans [Q (see also [|]) contains a corresponding result 
for moments of traces of random symplectic and orthogonal]] matrices (which they 
deduce using the work of Ram Q on Brauer algebras), which can be used to prove 
our theorems in half the range, that is the m-th moment of Zf is Gaussian if supp / 
lies in the interval (— l/2m, l/2m). We wish to have the full range so as to compare 
with zeros of quadratic L-functions, where linear statistics show mock-Gaussian 
behaviour in the same full range (this can be deduced from the work of Rubinstein, 
0). The case of Dirichlet L-functions, which correspond to the unitary group, was 
considered in ||. 

To obtain the results we desire, we abandon moments and instead use the cumulants 
Cf {N \g) oiTrg(U). These are defined via the expansion 

oo t 

logE GW (e^^)) = ^CfW( 5 )| 
l=x 

where Eg(jv) denotes the expectation with respect to Haar measure over the group 
G(N). The cumulants have previously been considered in this context by Soshnikov 
|s) (interestingly, his results again only give half the required range), and it is his 
combinatorial approach that we adopt. 

There is a natural decomposition for the cumulants on the symplectic and special 
orthogonal groups. For brevity we will describe the situation for the symplectic 
group (so TV, the matrix size, is assumed to be even). The cumulants can be 
written as 

C* p(N \g) = 2 l CtJU(g) ~ ^n+M 

We show that the odd parts C^ +1 (g) of the cumulants vanish in a certain region, 
and in fact if gk = for > (N + l)/£ then the £-th cumulant vanishes. 

For all <?, the even summand equals a unitary cumulant: 

CiT +1 {g) = \c^ N+1 \g) 

We may now employ the available results about the unitary group to deduce that 
Cf n+ a ^ so vanishes in a larger region. Setting g = we obtain Theorem ||. 



1 Notc that Diaconis and Evans consider orthogonal matrices, whereas we are interested in the 
special orthogonal group 
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Since moments and cumulants give essentially equivalent information, we can now 
go back to computing averages of the product of traces on classical groups and 
resolve a problem raised in jl], Remark 8.2], to show 

Theorem 3. Let Zj be independent standard normal random variables, and let 

1 if j is even 
if j is odd 



i) If aj G {0, 1,2,...} for j = 1,2,... are such that ^ jaj < N + 1, where N is 
even, then 

ii) If aj 6 {0, 1,2,...} for j = 1, 2, . . . are such that ^ jaj < N — 1 then 

E so(N) {n( Trt/3 ) aj } = E {l[(V]Zi + Vs)° i } 

Similar Theorems have been proven by Diaconis and Evans jl]] , though only for half 
the range (that is, they require ^2jaj < N/2). 

2. Cumulants of linear statistics 

In order to calculate C^ 9 ^ (g) we need to know the moment generating function. 
Weyl Q showed that E Sp ( Ar ){e* Tr£ '( c/ )} could be written as an integral over the 
N/2 independent eigenphases (recall that N must be even for a symplectic matrix 
to exist). He showed that, writing N — 2M, 

Esp(A0{e tTrs(L °} = Esp(JV) jexp 5 (0„)j j 

. M 

= / Dct{Q s P( 2M )(^,^)} 1 < l . J <MTT 

■AM" n=l 

where the kernel is Q Sp ( N '(x, y) := Sn+i{x — y) — Sn+i(x + y) with 

1 sin(iVz/2) 



e 2t ^ d9 n 



(1) S N (z) :-- 



2ir sin(z/2) 



Now, it is a general fact that if 6 n G T, where T is some real interval, are such that 

E | ex P (jy^n)] | = Det{g M (^,^)} 1 <i J < M n et9(6n) d0 « 
then if the £th cumulant of ^2g(0 n ), C'e, is defined by the expansion 
logE jexp (tJ2g(e n ^j j = £ t Cl 

then g § 

£ f rn 

(Xj)Q M (Xj,Xj +1 ) dXj 



4L. p lib 

( 2 ) ^ = E E (-ir + v-i)i / n^( a 

m=l a£P(l,m) Jm 1=1 
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where we identify x m +i with x±. Here P(£, m) is the set of all partitions of £ objects 
into m non-empty blocks, where the jth block has Xj = Aj(tj) elements (that is 
\ r .= #{i : l<i<l, ff(i)=J». 

Thus, 

C^ N \g) = 



£ „ rn 

m=l o-£P(.«,m) "'[0^] m J = 1 



Since Q Sp ^(x,y) is odd in both variables, Iljli Q Sp ^ N \%j, x j+\) is even in all 
variables, and so, since g is an even function, we may extend the integral to be over 
[—it, 7r] and thus 



C^ N \g) = 2^ E M) m+ V-l)!x 

m=l aeP(i,m) 
If m 

and on expanding out the middle product on the bottom line, 

c? w Cff)=2<x; E (-i) m+ v-i)^ E >< 

m=lrT£P(e,m) ei=±l,...,e m =±l 

J|.9 Aj ( a; i) e j'S'7V+i(^ - 

t -7r,7r] m 3 - =1 

= 2*Clffi_ 1 ( S )-2?C$ +1 ( S ) 

where Cf™ +1 {g) contains those terms with JXJLi e j = + 1 anc ^ ^Tlv+iG?) contains 
those terms with IIJli e j = — !• 

Similarly one can calculate the other groups, using Weyl's calculation of Haar mea- 
sure, which is summarised in table [l| 

2.1. Summary. Put 

£ 1 

(3) cs n (.9) = E E (- 1 ) m+1 ( TO - 1 ) ! ^r E 

m=la£P(l,m) ei=±l,...,e m =±l 

n<y=+i 



1J.'/ A :.S a;; •'•„ - ej-Xj+i) da 
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Group 


Tr g{U) 


Kernel Qm(x,u) 


Range T 


U(iV) 




Sn{x,v) 




(-7T,7r] 


Sp(AT) 
A 2.1/ 


2Ef=x5(^n) 


Sn+i(x -y) - S N+ i{x - 


\-y) 


[O.Tf] 


SO(iV) 
AT = 2M 




S N -i{x -y) + S N -i(x - 


r-tf) 


[0,7T] 


SO(iV) 
TV = 2M + 1 


s(O) + 2£f =lff (0„) 


Sn-i(x - y) - Sn-i(x - 


\-v) 


[O.Tt] 



Table 1 . Kernels for Haar measure over the classical compact groups 



and 



m=l a£P(e,m) 



(4) q#( ff ) = x; Yl (-ir+V-i)^ E 

ei=±l,...,e m =±l 
t -7r,7r] m =1 



with defined in ([!]). 
• For all t, 



• For all 



Cf (2M) ( 3 ) = - 2<C$l /+1 ( ff ) 



Q SO(2M) ( 5 ) = 2^is_ 1 ( 5 ) + s'q^.^j) 



• For £ = 1, 



CT {2M+1) {9)= 20^(3) -2Clf M {g)+ £ 5fc 

fc— — oc 

and for all £ > 2, 

c SO(2M + l) (5) = 2 ^cven _ 2 ^odd /(ff) 

In the next section, we will show that Cf^j(g) = hCf^ M \g), and then we will 
calculate Cg'j^(g), first in the case when M is odd, and then in the case when M 
is even. 
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The results will show that 

(5) C^ N \ 9 )=^^ N Hk 1 ,...,k l )l[9 kj 

k£Z e j=l 

where (i^ N \ki, . . . , kg) is invariant under permutations of its arguments. 
Combining the results from the next section proves the following theorems: 
Theorem 4. 

M 

-,Sp(2M) , 



oc oo M oc 

n=l fc=M+l 1=1 k=M+l 

and for i > 3, ftf^ih, ...,ki) = i/ELi \kj\<N+l. 



Theorem 5. When averaged over the special orthogonal group, the mean of Tr g(U) 
is 

M-l 

C 1 SO(2M) ( 5 ) = 2M 30 + 2^ ff2 „ 

n=l 

M oo 

C 1 SO(2M+1) ( 5 ) = (2A/+l). 90 + 2^.g 2 , i + 2 £ <?„ 

n=l n=2M+l 

and £/ie variance is 

oc oo M-l co 

-,SO(2M), 



C 2 so(2M) («,) = 4 £ mm(n, 2M - l)g 2 n + 4 £ 5 ^ + 8 £ £ 

n=l fc=Af (=1 fc=Af 

oo 2M-1 oo 

Cf ° (2M+1) (9) = 4 min(n, 2M)fl£ - 8 9{m+n)/l9{m-n)/2 



n=l n=l m=2M+l 

n odd m odd 

For£>3, t i S e O(N) (k u ...,k e )=0 i/E^il^l <N-1. 

3. The combinatorial calculations 



3.1. The calculation of Cf^{g). The following lemma was stated by Soshnikov 
in §: 

Lemma 6. For all I, 



Ct:^{9) = \c u t [M) {g) 



Proof. Symbolically denote 



/lib 
'Y[g x >(x 3 )S M (3 



(6) / i f g Xl (xj)S M (xj - ejXj+i) dx 



3 C J^J + -W ^3 

3 = ' 



by (ei, 62, ... , e m ). If ei = 1 do nothing, but if ei = —1 then change variables to 
X2 1— > —2:2, and note that since g and S"m are even functions, and the integral over 
X2 is over [— 71", n], then (^) becomes (+1, —£2, £3, • ■ ■ , £m)- 
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Observe that this achieves the following: If the initial situation was (—1, —1, . . . ) 
then it becomes (+1, +1, . . . ) while if it was (— 1, +1, . . . ) it becomes (+1, — 1, . . . ). 
Therefore there is either the same number of — l's in the set of e or there are two 
less —l's. 

Now repeat for the new £2, changing variables only if it is —1, and so on all the 
way up to e rn . Each time the action either leaves the number of —l's unchanged or 
reduces it by 2. Since we started with an even number of —l's in the set of e this 
algorithm will terminate with (Q) equaling (+1, +1, . . . , +1), which is independent 
of e. There are 2 m ~ 1 possible e with an even number of —l's, and so 



Cf,™ = £ E (-i) m+ V-i)!^r>< 

m=l aeP(l,m) 

\{g^{x )S M {x -x j+1 ) dr, 

which we recognise as ^C^ M \g). □ 

The cumulants of a random unitary matrix have previously been calculated, essen- 
tially by Soshnikov J8|, but they can also be deduced from the work of Diaconis 
and Shahshahani [|| and of Diaconis and Evans . 

Theorem 7. (Soshnikov). Let be the £th cumulant of Tr g(U), averaged 

over all N x N unitary matrices with Haar measure. Then 

= Ngo 

OO 

c u(N) _ ^2 min(|n|, N)g n g_ n 

and for £ > 3, 



n— — 00 



C^ [N \g) < consk ^ l fc illfffcil ■ ■ ■ Iff** I 

feiH \-k e =0 

|fei|+-+|fe«|>2iV 

Remark. The heart of the proof of this theorem is a deep combinatorial fact called 
the Hunt-Dyson formula. 

Remark. Actually, the error term in Q has the sum running over all fci + - • - + kt — 
such that \k\\ + ■ ■ ■ + \k^\ > N . But it is clear from equation 2.9 of || that there is 
no contribution to for £ > 3 if ^ ki1^ i>0 y < N and if ^ — kil^ ki<Q y < N. 

Since the ki sum to zero, it must be that the sum over positive terms equals the 
sum over negative terms, and so this is the same as the condition that \kj\ < 2N, 
as we have it in the theorem. 

3.2. The calculation of C^ d M+1 (g). Observe from (Q) that 



(7) S 2M +i{z) = — E 



M 

I 'l l z 



e 

2vr * 
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Lemma 8. One can calculate C°^ M+1 (g) and 2M+i(s) exac ^V- 

1 M 

C°,2M+i(.9) = 2 E 52n 
n=-M 

1 M 

Q>,2M+i(-9) = 2 E X! 9l+k9l-k 

l=-M \k\>M 

Proof. First of all, from (Q) we have that 

1 r 



C^ M+1 (g) = ^ / g(x)S 2M+ i(2x) dx 



C° 2 f M+1 {9) = \T g 2 (x)S 2M+1 (2x) dx 



g(x)g(y)2S 2M +i(x + y)S 2M +i(x - v) dx dy 



and using (Q) we see that 

1 M 

Cl,2M+l(.9) = ; ^ 92r. 



n=-M 

and 

A/ oo , M M 



1 . 1 

c 2 od 2 d M +1 (. 9 ) = ^ E E 

l=-M k=-oo l=-Mk=-M 
M 

= ~ E E 



gi+kgi-k 



2 

Z=-M |fc|>M 

as required. □ 
Lemma 9. For £ > 2, 

| C°,2M+i (5) | < eoras^ E Ife I • • • Ifffc* I 

kez £ 

|fci|H h|fe«|>2M+l 

Proof. Fix er 6 P(i, m), and for k = (ki, . . . , kg) 6 1> e set 

Ai 

A'i = ^> 

1=1 
A1+A2 

J=Ai+l 



(=A 1 +---+A m _ 1 + 1 
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(recall that I = Ai H h A m ). Therefore 



1I'/ X "V XII'/ II' K 



Hence, the integral in (^) 

JJff A ^(a: 3 -)S2M+i(a;j - e j x j+i) dx j 

t -7r,7r] m 3 - =1 



E EII^/ J^e i ^^e i ^^- e w)^ 



2- 



-M<m,..,n m <Mk6Z'!=l ■> [-' tt.ttJ J = 1 

2vr 



/in 
n ° xp (^(^ + n 3 - e i-i n J-i)) 
kgi- 1 = 1 -iHsni,...,n I „sM L- 7r ^] m j = i 

where we have used (]?]) to express S^M+iC^j ~ e j x j+i) hi hs Fourier representation, 
and we have defined £q = e rn , no = n m (so all indices are cyclic). 

The integral above will be 1 or depending on whether rij — ej-irij-i = —Kj or 
not, so defining 

(8) Af{M, a, k, e) = 

# {-A/ < ni, . . . ,n m < M : rij - ej-Wj-i = -Kj , j = 1, . . . ,m} 
(the i^i, . . . , -Km, depend on both k and a, recall) we see that 

(9) Q^ + i(p)=EII^E E ( ~ ir+ 2 1 ,i m " E -^(M,a,k, £ ) 

keZ*J=l TO=l<reP(f,m) ei,...,e m =±l 

n^=-i 

Lemma 10. Let YYjLi e j = ~ !• ^ e?1 M(M, a, k, e) is either or 1. 



If ^2 1— i k\ is odd then Af(M, o~, k, e) =0. 

If^2i—i ki is even and ^/=i 1^(1 — then N{M,o~,\t,e) = 1. 



(proof deferred until the end of this section) . 
Therefore, if ^Li l fc * < 2M + 1 then 

( 10 ) E E (- i r +1 ( m ~ 1 ) ! ^r E ^(M,a,k,e) 

m=l o-EP(i,m) ei,...,e m =±l 

n <y=-i 

= ^(k)E E (-l) m+1 (™-l)! 

m=l aeP{t,m) 

where 



X(k) = 



f if ^/ * s even 

otherwise 
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Using the fact that for I > 2, 

E E (-l) m+1 (m-l)! = 

m=l aeP(£,m) 

we see that (jic| ) vanishes for Yli=i — %M + 1 if £ > 2, Inserting this into (||) 
and estimating the contribution from the terms with Ym=i ^ 2M + 2 we see 
that 



|Q° d 2M + i(.9)| < const, E IT 1^ I 

£f =1 |fc,|>2M+2 

This completes the proof of Lemma |9[ □ 



Proof of Lemma 
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We treat all indices as cyclic modulo m. So no — n m and n m+ \ = n\ etc. 

We assume that Yl^Li e j = — 

Define the to x to matrix to be such that 



6i_i if j = i — 1 
otherwise 



so that 

(#n)j = ej-mj-i 

From the definition of Af(M, a, k, e) (which is given in (||)) we see that it is the 
number of solutions of (I — E)n = — K subject to —M < nj < M. 

Now, 

= e J-l £ i-2 ■ ■ • £j-k n j-k 

and so E m — e\ . . . e m I = — I by cyclicity of indices and the assumption that 

Hence 21 = 1 - E m . But 7 - factorizes as 

I - E m = (I - E)(I + E + ■ ■ ■ + E rn - 2 + E m - 1 ) 

and therefore 

(I - E)- 1 = ~(I + E + --- + E m - 2 + E" 1 - 1 ) 

If we ignore the restriction that — M < rij < M then, over the reals, there is exactly 
one solution to (I — E)n = — K which is 

(11) 

Tlj = — -z (Kj + Cj-iKj-i + ej-i£j- 2 Kj-2 + • • • + Cj-iCj-2 ■ ■ ■ £j-m+lKj- m +l) 
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This is a solution over the integers if rij is an integer, which will be the case when 
the term inside the bracket is even. Since £j = 1( mod 2) for all j, the term inside 
the bracket is even when 

m t 

Kj + Kj-i + ■■■ + K j - m+1 = J^ifc = ^ *j 

»=i i—i 

is even. There are no solutions over the integers when this is odd. (Note that the 
even and oddness is independent of e and of the partition a). 

Finally, one must check that the condition — M < rij < M holds. From (|ll]) we see 
that 

K|< -El^l< 2^ |fci| 

i=l i=l 

and so if we assume that Yjt=i l^'l — ^Af, then the condition holds. 

Thus Af(M, a, k, e) =0 if 2;=i i s °dd, and M(M, a, k, e) = 1 if Y^i=i is even 
and Yfi=i \H < 2M. 

This proves Lemma |l0[ □ 



3.3. The calculation of (.<?)• Basically, this section is like the previous, with 

the essential change being that 

2M-1 

S 2M {z) = ^ £ e-^ 

n=-(2M-l) 
n odd 

as opposed to (]?]) which says 

1 2M 

S 2M+1 (z) = - £ e-^/ 2 

Z7T ' — ' 

n=-2M 



Lemma 11. One can calculate C°^f M (g) and C% 2m(9) exa ctly. 

1 M 

C°,2M(9) =2 X] ff2n-l 
n=-(M-l) 



2A/-1 



^2,2m(s) — 9 £ £ 9 , A(n+m)fl'i(n-m) 



:-(2M-l) |m|>2M+l 
n oc£(Z m odd 



Lemma 12. For £ > 2, 



| C°,2M (5) I < const e £ |^ fcl | . . . \g k 

|fei|H h|M>2M 
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The proof goes through the same as before, with equation (^) becoming 

Nodd{M, a, k, e) = #{-(2M - 1) < rij < 2M - 1 , rij odd : 

\n 3 - ej_i|nj_i = -Kj , j = l,...,m} 
Rewriting equation (0) we see the solution requested by N dd(M, a, k, e) is 

rij = - (Kj + e-j-iKj-i + ( , : < , 2 l\ , 2 + • • • + ej_i£j_ 2 . . . ej-m+ii^-jn+i) 

so long as rij is odd and — (2M — 1) < rij < 2M — 1 (and there is no solution 
otherwise) . Therefore Lemma [lO] becomes 

Lemma 13. Let YYjLi e j = ~ !• Nodd{M, a, k, e) is either or 1. 



If ^2 1— i ki is even then N dd{M, k, e) = 0. 

// Y2i=i ^1 is odd and Ym=i ^ I — — 1 i/ien M dd{M, a, k, e) = 1. 



4. Moments of traces 

We will now use Theorem ^| to prove the second part of Theorem ^. (The proof of 
the first part from Theorem [| being analogous). 

Recall from (0) that 

e 

c $o(n) ^ _ fJ s i 0( - N) (m , . . . , m) J| g nj 

n£Z e 3 = 1 

where ii^ ^ N \ni, . . . , ng) is invariant under permutations of its arguments. Assum- 
ing go — 0, then we have 



• If |ni| < N then 



^SO(JV), n _ I 1 if m ^ is even 
1 | otherwise 



• If |m| + |n 2 | < N then 

SO(JV), \ 
M 2 («4,™ 2 ) 



|m| if |m| = \n>2\ 

otherwise 



• If £> 3 and £*. =1 |rij| < N then ( J (m> . . . , n e ) = 0. 
It is also true that if go = 0, 

E G {(Tr g(U)-C?(g)) m } 

m 

(12) =2 m ^ E G {(TTU^-^(n 1 ))...(TrU n --^(n m ))}Y[g nj 

n£N m j=l 

( n) _^(cf(9)\ k2 (cng)\ k3 (c? n (g)\ km mi 

1 1 ^\ 2\ ) \ 3\ ) ml J fc 2 !fc 3 !...fc m ! 

where the second sum runs over all values of kj > such that -_ 2 j'fcj = m (it is 
simply writing the mth moment in terms of its cumulants, having subtracted the 
mean) . 
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Let a,j £ {0,1,2...} for j = 1, 2, . . . by such that ^ jaj < N. Define 

Vj = 



1 for even j 
for odd j 



so that (J*i° (J) = r/j for |j| < N. 

Putting m — Y^, a ji we will evaluate the coefficient of Ilfe)^ m ( P^ ) an d m 
the two being equal to each other. 

Consider first equation (|T^). The coefficient of n(5j) aj m 

m 

2 m J2 Eso(iv) {(Tr C/" 1 - r? ni ) . . . (Tr U n ™ - Vn J} [] g nj 

equals 

(14) __E SOW {n(Tr^-^r} 

Consider next equation (|l3| ) . Note that the restriction on the aj means that there is 
no contribution to the coefficient of nfe)^ fr° m C^°^ (g) for all £ > 3. Therefore 
the coefficient in (Ol) is if m is odd and is the coefficient of Il(ffi) ai m 



,„,3 • 



2 m /2( TO /2)! V 2 W V 2™/ 2 (m/2)! 
if m is even. This coefficient is zero unless all the aj are even, in which case it is 

(151 m! 2 m/2 (m/2)! TT7 a ' /2 

( } (m/2)! 2 n(a,72)! 11J 

(to see this, note that the structure of /i^ ^ means that n-ij must equal n-ij-x for 
j = 1, . . . , m/2. The second prefactor is just the number of ways of picking m/2 
integers such that ai/2 of them equal 1, 02/2 of them equal 2 etc.). 



Setting (|14|)=(15) and recalling that m = Y a j > we have 

mi f Ff i a ^ 2 a./2°/^ ! ; s, if all the a,- are even 
(Tr^- % ) a 4 = ^ U 2^/ 2 (a J /2)! i 



otherwise 



where Zj are iid normal random variables with mean and variance 1. 
Observe that this can all be rewritten as 

and is valid so long as X)i a i < N. 
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